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Abstract 

Spin relaxation of conduction electrons in metals is significantly influenced 
by the Fermi surface topology. Electrons near Brillouin zone boundaries, 
special symmetry points, or accidental degeneracy lines have spin flip rates 
much higher than an average electron. A realistic calculation and analytical 
estimates show that these regions dominate the spin relaxation, explaining 
why polyvalent metals have much higher spin relaxation rates (up to three 
orders of magnitude) than similar monovalent metals. This suggests that 
spin relaxation in metals can be tailored by band-structure modifications like 
doping, alloying, reducing the dimensionality, etc. 
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The ability of normal metals to carry a spin-polarized current has led to the develop- 
ment of spin transistor Q, a ferromagnet-normal metal-ferromagnet sandwich device that 
can switch current depending on the relative orientation of the magnets. Spin polarized 
currents also flow in the metallic layers of giant magnetoresistance heterostructures || , sub- 
stances promising for magnetic signal processing. A major limitation of the quality of these 
devices is imposed by the time (the so called spin relaxation time T\ ||) an unbalanced spin 
population in metals persists; the metallic layers used in these devices cannot be thicker 
than the spin diffusion length which is proportional to \fT[. Finding ways of enhancing Ti 
is technologically important also from the perspective of quantum computing that represents 
bits by electronic spins [[|. One recent successful attempt to increase Ti in semiconductors 
by doping is described in Ref. [f|. 

We recently reported || on a theoretical study of the spin relaxation in polyvalent metals 
where we showed how the band structure affects T\. The physical picture is the one of a 
random walk on the Fermi surface: the weakness of the spin-flipping interactions in a typical 
metal ensures that an electron changes its momentum many times (typically a ten thousand) 
before its spin flips. At some points (we call them spin hot spots), however, the spin-flip 
probability is much enhanced. Near a Brillouin zone boundary (Bragg plane), for example, 
this probability can increase a hundred times. If the electron jumps in or out of a special 
symmetry or accidental degeneracy point, the probability may be even close to one, that is, 
the spin-flipping and spin- conserving jumps may be almost equally frequent! Although an 
electron jumping into a spin hot spot is a rare event, it may nevertheless dominate the spin 
relaxation. 

The above picture solves a longstanding experimental puzzle formulated first by Monod 
and Beuneu 0: Why spin in some metals decay unexpectedly fast? Experiments show that 
alkali and noble metals have spin relaxation times consistent with the predictions of the 
theory of Elliott || and Yafet f9j, when the spin-flip rates are calculated from the atomic 
state parameters. Metals Al, Pd, Mg, and Be, on the other hand, have spin relaxation 
times much smaller (up to three orders of magnitude for Mg and Be) than estimated. This 
discrepancy is disturbing since it shows that otherwise similar metals may have very different 
spin relaxation times. Silsbee and Beuneu [T(| were first to notice that accidental degeneracy 
points can significantly increase the spin relaxation of aluminum. We substantiated this idea 
by a rigorous calculation, and extended it to include all band-structure degeneracies. The 
experimental puzzle is then resolved by recognizing the fact that the two groups of metals 
have different valency. Alkali and noble metals are monovalent-their Fermi surface lies 
completely within the first Brillouin zone. The spin-flip rates are more or less uniform with 
the values close to the ones derived from atomic physics (there are no spin hot spots). On 
the other hand polyvalent Al, Pd, Mg, and Be have numerous Fermi surface anomalies 
were spin hot spots can be formed. Aluminum, for example, have spin relaxation rates 
determined primarily by the Fermi surface regions near the Brillouin zone boundaries and 
around the accidental degeneracy points P,|T0| . In beryllium, the Fermi surface cuts through 
a degeneracy plane! 

Since no polyvalent metals other than Al, Pd, Mg, and Be have been measured for 
Ti so far, our explanation of the experimental puzzle becomes a prediction for future T\ 
measurements. In addition, the spin-hot-spots picture suggests a way of altering Ti via 
band-structure modifications. Doping into the conduction band, for example, may shift the 
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Fermi surface away from some special symmetry points and increase T\. Similarly, alloys 
or systems with reduced dimensionality will have spin relaxation rates different from those 
of the corresponding elemental or bulk metals, respectively. A general rule of thumb for 
increasing T\ of electrons in a conduction band is washing out the spin hot spots from the 
Fermi surface. In what follows we introduce the basic concepts of our theory and discuss 
the results in a qualitative fashion. 

If the periodic potential due to ions in a crystal lattice contains spin-orbit coupling (a 
term proportional to the scalar product of the orbital and spin momentum operators, L- S), 
the electronic Bloch states are a mixture of spin up |f) and down ||) species ||: 

*L(r) = kn(r) |T) + 6 kn (r) exp(*k ■ r), (1) 
= K kn (r) ||) - b\ n (r) |T>] exp(*k ■ r). (2) 

The lattice momentum and band index are k and n, respectively, and dk n (i") and &kn( r ) 
are complex periodic functions with the period of the lattice: if G denote the reciprocal 
lattice vectors, then ak n (r) = J2g a k«(G) exp(iG • r) and similarly for bk n (r). Both states 
have the same energy E^ n , as follows from time and space inversion symmetry ||; the 
numbering of bands is therefore the same as without the spin notation. The degenerate 
states \I/ kn and \l/ kn are chosen to represent electrons with spins polarized along z direction 
@: (^kJ^l^jL) — — (^icnl^l^kn) < and the off-diagonal matrix elements are zero. 
This condition implies that a^ n (r) have values close to one, while &kn(r) are much smaller, 
decreasing with the decrease of the strength of the spin-orbit interaction (with the exception 
of the points where the spin-orbit interaction lifts a degeneracy). 

Elliott P| pointed out that ordinary (spin conserving) impurity or phonon scattering can 
induce transitions between \l/ kn and ^l_, n , (either their spin up or down amplitudes), leading 
to the flip of a spin polarization and thus spin relaxation. Simplifying assumptions lead to 
the formula f| 

1/H « 4(fo 2 )/r, (3) 
where (b 2 ) is the Fermi surface average of 

|&kn| 2 = £lMG)| 2 , (4) 
G 

and r is the momentum relaxation time. The weakness of the spin-orbit interaction makes 
the average spin-mixing parameter (b 2 ) much smaller than one and 1/Ti <C 1/r. At low 



temperatures, in a very pure sodium, for example, T\ can reach a microsecond ||11|| , much 
larger than momentum relaxation times which would reach a fraction of a nanosecond in 
similar samples. Note that Eq. ^| implies a series "spin resistor" model (scattering rates 
are additive) for the spin relaxation, if |&kn| 2 / r i s interpreted as the spin-flip rate for the 
scattering from, or to state kn. The reason is that because 1/71 ^ l/ r , a single electron 
experiences many different ("series") scattering events before changing its spin. Therefore 
the Fermi surface states with the highest spin-flip rates count most. (A counterexample is 
conductivity which is monopolized by the states with the lowest scattering rates.) 

In addition to the Elliott's mechanism, impurities and phonons affect the spin relaxation 
in other ways. Impurities induce the spin-orbit interaction which allows a direct transition 
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between the spin up part of \& kn and the spin down part of x I / j c / n /. The resulting spin 
relaxation is independent of temperature (as is the Elliot's impurity relaxation) and can 
be experimentally controlled. Similar transitions can be induced by the phonon-modulated 
spin-orbit interaction. This intrinsic effect leads, in principle, to a spin relaxation that is 
as effective as the Elliott's phonon- induced spin relaxation, but as we show below, becomes 
unimportant in many polyvalent metals where the Elliott's mechanism is enhanced by band 
structure. 

In most cases the spin-orbit interaction can be incorporated into the band structure 
as a perturbation, leading to the estimate (b 2 ) ~ (A/A) 2 . Parameter A is some effective 
spin-orbit coupling and A is a typical (unperturbed) energy difference between neighboring 
bands. Monod and Beuneu J7| estimated both A and A from their atomic values and found 
widely varying results for different metals. While this atomic substitution works well for 
alkali and noble metals, other metals (Al, Pd, Mg, and Be) seem to have (b 2 ) much larger 
than calculated. The case of aluminum and sodium is instructive. These metals have 
similar atomic numbers so their atomic A/A are also similar (within 10% |?J). Yet their 
spin relaxation times differ by two orders of magnitude P,p^-[T5|! The reason is the very 
different band structure of the two metals that gives very different values for A. In sodium, 
A is of order Ep, the Fermi energy. In aluminum, A varies between zero (here perturbation 
theory for degenerate states gives \b^ n \ ~ 0.5) and E F . Remarkably, the regions with 
small A (the spin hot spots) have enough weight to significantly increase the average (b 2 ) 
beyond the naive estimate (X/Ep) 2 which works so well for monovalent metals. A solid state 
environment greatly influences (b 2 ) and Ti, leading to what we call the band renormalization 
of the spin-orbit mixing (b 2 ). 

Figure [l] shows our calculation |J of the distribution p of the values attainable by |&kn| 2 
over the Fermi surface of aluminum. The width of the distribution is impressive-almost 
seven decades! Most states have |&k«| 2 below 10~ 5 . Higher values are much less frequent 
but they stretch up to about 10 _1 (we lacked enough precision to reach the upper limit 
0.5) with p decreasing linearly, as seen in the inset. This long tail, however, has a marked 
influence on the average (b 2 ) ~ 2.0 x 10~ 5 which is an order of magnitude larger than the 
value where p is maximal. The unusual character of the distribution is further emphasised 
by the plot (also in Fig.|Tj) of p for a hypothetical case of monovalent aluminum. The band 
structure parameters (like the lattice structure or the form of the electron-ion potential) 
are unchanged so that any difference in p between the two metals is caused solely by their 
different valencies, and therefore different Fermi surface geometries. The distribution of b 2 
is now relatively narrow with the average (b 2 ) ~ 3.4 x 10~ 7 coinciding with the center of the 
distribution; this value is about fifty times smaller than the average for (trivalent) aluminum. 
Adjusting for the density of states (monovalent aluminum would have 1/r reduced ~ 3 1 / 3 
times) the spin relaxation would be about seventy times slower. This is about the measured 
difference between sodium, which has similar Fermi surface as monovalent aluminum, and 
aluminum ||. 

The difference between the two very different forms of p in Fig. [I] is clearly attributed to 
the subtleties of the band structure of aluminum. The Fermi surface of monovalent aluminum 
is nearly spherical and lies entirely within the first Brillouin zone; the protuberances near 
the zone boundaries do not touch the planes as in the case of the noble metals. On the other 
hand, the Fermi surface of aluminum crosses the Brillouin zone boundaries and accidental 
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degeneracy lines [HJ. The distorted electronic states near the crossings are responsible for the 



long tail of p. Consider a band structure computed without the spin-orbit interaction. The 
spin-orbit interaction, being a part of the periodic lattice potential can induce transitions 
only between states whose k differ by a reciprocal lattice vector, or, in a reduced-zone 
scheme, between the states with the same k but different n. Let, for a state kn on the Fermi 
surface, the closest band to n is separated from n by A. The spin-orbit interaction mixes the 
spin amplitudes from the two bands (more distant bands to n will have smaller contribution 
and we neglect them) leading to the spin-mixing parameter |&kn| 2 ~ (1 — A/VA 2 + 4A 2 )/2. 
Three cases can occur. (A) For a general point on the Fermi surface, the band separation 
is of order E F , the Fermi energy, so that A ^> Vgo and |&k«| 2 ~ {Vsol Ef) 2 ■ This is the 
case of monovalent aluminum: parameters A ~ 3 meV (same as for trivalent aluminum ||) 
and Ep pa 6 meV give |&kn| 2 ~3x 10~ 7 (~ 10~ 66 ) in accordance with Fig. [1|. In (trivalent) 
aluminum typical spacings between bands are somewhat smaller than Ep ~ 12 eV; the 
estimate of A pa 3 eV is in reasonable agreement with Fig. [1]. (B) If the state is close to a 
Brillouin zone boundary associated with G, the band separation is ~ 2Vg (Vg is the Gth 
Fourier coefficient of the non-spin part of the lattice potential). Since typically Vg ^> Vso, 
|&kn| 2 ~ (Vso/2Vg) 2 ; this can be a few orders larger than in (A). Aluminum has V m ~ 0.2 
eV which gives |&kn| 2 about 6 x 1CT 5 (~ 1CT 4 ' 2 ) coinciding with the onset of the tail in Fig. 
[l]. Curiously, the Fermi surfaces of the noble metals too come into contact with some zone 



boundaries. The noble metals, however, have unusually large Vg []17] so the estimate (A) 
works equally well. Finally, (C) the spin-orbit interaction can lift the degeneracy of two or 
more bands. The mixing of spins is complete and |&kn| 2 ~ |«kn| 2 ~ 0.5. Even this case occurs 
in aluminum, where the second and third band accidentally coincide at some Fermi surface 
points (the accidental degeneracy points). The neighborhood of these points is responsible 
for the long tail of p in Fig.|I|. The spin hot spots are the states with the properties (B) and 
(C). 

As Fig.[TJ reveals, the spin hot spots essentially determine the average (b 2 ). We give a 
simple example demonstrating that this is reasonable. Consider electrons in the presence of a 
single Brillouin zone plane that is associated with the reciprocal vector G. If the electron- ion 
potential (or, rather, pseudopotential) is weak, the band structure of such a system is found 
by using just two plane waves, exp(ik • r) and exp[i(k — G) • r]; this is sometimes known as 
the two orthogonalized-plane-wave (OPW) method [TjJ. The problem has axial symmetry 



along the direction of G, so we conveniently shift the origin of k to G/2 and decompose k 
into the components parallel (z) and perpendicular (r) to G: k = G/2 + r + z. To further 
simplify the notation, energy and momentum will have the units of h 2 (G/2) 2 / (2m) and G/2, 
respectively; m will be the electron mass. The band structure is formed by two energy levels 



with the dispersion (see, for example, Ref. [18]) 
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E ± = l + r 2 + z 2 ±^sgn{z)A{z), (5) 



A(z) = ^IQz 2 + 4V<|. (6) 

We allow z which measures the distance from the plane, both negative and positive values 
(repeated-zone scheme) and take E + as the reference level. The Fermi surface is then the 
revolution of the curve 

E + (r,z) = E F (7) 
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about G. If z is negative (states inside the first Brillouin zone) E~ is the upper band, if z 
is positive E~ is smaller than E + ; the spacing between the two bands A(z) depends only 
on z. 

The spin-orbit interaction will mix the spin amplitudes of the two bands. In most cases 
A < Vg and we can use the estimate (B): 

K n \ 2 =\b z \ 2 ^[\/A(z)] 2 . (8) 

A proper evaluation of |&kn| 2 would, instead of effective A, have the matrix elements of 
the spin-orbit interaction between the band states with energies E~ and E + . Such matrix 
elements, in general, depend on both r and z, and for spin quantized not along G, also 
on the angle of revolution <fi. In particular, the spin mixing would vanish at k = G/2 
(that is, r — z — 0). This follows from group-theoretical arguments or, more simply, 
from the fact that the spin- mixing matrix elements between the plane waves exp(zk ■ r) and 
exp[i(k — G) • r], out of which the band states are formed, are proportional to the vector 
product k x (k — G); this product vanishes for k = G/2. Incidently, this vanishing of spin 
mixing is another reason why the band structure does not affect T\ for the noble metals 
despite the fact that some of their states come into contact with zone boundaries at G/2. 
Introducing effective A therefore overestimates the region around G/2. In most interesting 
cases of polyvalent metals, however, the Fermi surface crosses the plane far enough from this 
point (that is, if z — radius r is a significant fraction of G/2) and our approach is justified. 
The Fermi surface average of | b z | 2 given by Eq. |8] is 

1 rz2 A 2 

The limits of the integration, z\ and Z2, are the minimum and maximum value of z reached 
by the Fermi surface; they are obtained by putting r = and solving Eq. |7| for z. We 
are interested in how (b 2 ) depends on the Fermi surface geometry. The Fermi energy which 
determines z\ and z 2 will now be a variable parameter (simulating doping). Although 
the integral in Eq. |9] can be evaluated analytically it is more instructive to consider some 
limiting cases only. The numerical estimates will be done with the band-structure parameters 
appropriate for aluminum so that the results can be directly compared with the full band 
structure calculation of Fig.|I|. The parameters are A ~ 4.3 x 10 -4 and Vg = Vm ~ 0.035 
(we now use the dimensionless units for energy). Although the aluminum Fermi surface 
crosses both (111) and (200) zone boundaries, we choose as Vg the Fourier coefficient Vm 
which is much smaller than V200 and has therefore greater impact on (b 2 ). 

If Ep is small, the Fermi surface is a small, almost undistorted sphere. Both z\ and z 2 
are close to —1 so we can neglect Vg in Eq. || and obtain (b 2 ) « A 2 /16 ~ 1.2 x 10~ 8 . This 
number is not affected by the presence of the Brillouin zone boundaries and can be thought 
of as an atomic limit for the spin mixing. When Ep is large enough for the Fermi surface 
to cross the zone boundary, z 2 > — Vg and the small region of the states around the plane 
gives the dominant contribution to (b 2 ). Indeed, when taking the integral in Eq. ^| from 
— Vg to Vg only, A(z) w 2Vg and (b 2 ) w A 2 /(4Vg&f); the Fermi momentum k F = \fE~p~ 
is given by z 2 — Z\ ~ 2kp. Before substituting numerical values we note that there are 
eight (111) planes in aluminum. They contribute independently to (b 2 ) giving the estimate 
(b 2 ) ~ 8 x 10 -6 . The full calculation in Fig. [I] gives the value which is more than twice 
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larger. This is explained partly by our neglect of the (200) planes and, more important, by 
ignoring the accidental degeneracy points. 

Figure |2] plots the full dependence of (b 2 ) on Ep within the two OPW approximation for 
the eight (111) aluminum planes. Doping, or increasing the size of the Fermi surface can 
increase (b 2 ) by almost three orders of magnitude! As the Fermi surface expands towards 
the zone boundaries, the Bloch states that are closest to the planes (they have smallest 
A) become more and more important. Once the Fermi surface reaches the plane (at the 
energy 1 — Vq ~ 0.97) the average (b 2 ) goes through a maximum and then slightly decreases. 
The decrease is tied to the saturation of the integral in Eq. £| only the states in the small 
neighborhood of the plane are dominant and (b 2 ) ~ l/(z 2 — £i)- The quantity z 2 — zi is up to 
a numerical factor the density of states. As Fig. [l] shows, in this toy calculation monovalent 
aluminum would have (b 2 ) about fifteen times smaller than aluminum, in good agreement 
with the results of Fig. |]. 

The qualitative understanding of the above estimates for the impact of a Brillouin zone 
boundary on (b 2 ) is quite simple. We need to estimate the probability with which an electron 
in its random walk on the Fermi surface jumps into a state close to a Brillouin zone boundary. 
Alternatively, we can ask what is the portion of the free-electron Fermi surface states that 
have A < 2Vg (these states will be significantly perturbed by the presence of the boundary). 
The answer is ~ Vq and the reason that this number is linear in Vq (the linearity is crucial) 
is that A increases linearly away from the plane, as in Eq. |6|. The average value (b 2 ) is 
then ~ Vq X A 2 /^ ~ ^ 2 /Vg- The enhancement of the spin mixing due to the presence of 
a single Brillouin zone boundary is then ~ 1/Vq, as in Fig. 0. 

To estimate how accidental degeneracy points affect (b 2 ) is more difficult [|TD | . Consider 
two bands coming into contact at a single point R, different from the symmetry points, 
on the Fermi surface. As noted in Ref. |T3 the band spacing A grows linearly as we go 



away from R in almost all directions (the exception is the direction along the corresponding 
accidental degeneracy line). We need to divide the region around R in two: the first has 
states with ^ A < 2A, the second with 2A ^ A < 2Vg- Since A ~ 8k, where 8k measures 
the distance from R, there are ~ A 2 points in the first region, all with (approximately) the 
same value of &kn ~ 1 from estimate (C). The average (b 2 ) therefore scales as A 2 . The second 
region is different in that we have b^ n depending on A [estimate (B)]. The average (b 2 ) is 
thus proportional to 



(6 2 )^ /%AA^. (10) 



This integral evaluates to A 2 1u(Vg/A) giving a weak enhancement of about ln(Vc/A) ~ 4. 
More accurate evaluation based on a four OPW approximation would give an additional 
enhancement of (V^oo/^in) 2 ~ 10 [§[. As Vq becomes smaller the importance of accidental 
degeneracy points diminishes. 

Another type of a spin hot spot, not relevant for aluminum though, is the region around 
a special symmetry point. It may happen that the Fermi energy coincides with a set of 
two or more degenerate levels at a symmetry point. This is the case, for example, of the 
fee palladium and platinum, whose Fermi surfaces go through the fee L points | 20| . If the 



spin-orbit interaction lifts this degeneracy, the renormalization of (b 2 ) can be significant. 
This effect is, however, not easy to estimate qualitatively. For the fee W point we find that 
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the enhancement of (b 2 ) relative to the atomic value of ~ A 2 is about Vg/A ||19|| . This can 
range from a ten to a thousand. 

Our final note concerns the hexagonal Mg and Be, where the deviation of T% from the 
simple estimates is most striking We argue that this is also a manifestation of the band 
renormalization of (b 2 ). Without the spin-orbit interaction, all the states at the hexagonal 
faces of the first Brillouin zone of a simple hexagonal structure are degenerate [^TJ. The 
spin-orbit interaction lifts this degeneracy P| (except at some symmetry points and lines), 
presumably by the amount VgA, the largest second-order term containing the spin-orbit 
interaction (any first order term vanishes since the structure factor associated with the 
hexagonal faces is zero []2l|). The contribution to (b 2 ) of the points where the Fermi surface 
intersects the hexagonal faces is ~ VgA: the characteristic value |&kn| 2 ~ 1, times the area 
of the affected part of the Fermi surface, VgA. The enhancement measured in terms of A 2 
is then V"g/A; this can be as large as a thousand for light elements like Mg and Be. 

We thank P. B. Allen and M. Johnson for helpful discussions. This work was supported 
by the U.S. ONR. 
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FIG. 1. Calculated distribution p (in arbitrary units) of the spin- mixing parameters |6kn| 2 for 
aluminum. The corresponding average (b 2 ) ~ 2.0 x 10~ 5 is indicated by a solid arrow. The linear 
tail of the distribution is shown in the inset. The dashed line shows what the distribution would 
be if aluminum were monovalent ((b 2 ) « 3.4 x 1CP 7 , dashed arrow). 



FIG. 2. Two OPW calculation (with the parameters Vq and A suitable for aluminum) of the 
average spin-mixing parameter (b 2 ) as a function of the Fermi energy Ep. The values of (b 2 ) 
that correspond to monovalent (open circle, Ep = 0.81) and trivalent (filled circle, Ep = 1.7) 
aluminum are 4.1 x 10 -7 and 6.2 x 10 -6 , respectively. The vertical line indicates the case where 
the free-electron Fermi sphere comes into contact with the Brillouin zone boundaries. 
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E F (reduced units) 



